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Whether it occurs in superconductors, helium-3 or inside a neutron star, fermionic 
superfluidity requires pairing of fermions, particles with half-integer spin. For an 
equal mixture of two states of fermions ("spin up" and "spin down"), pairing can be 
complete and the entire system will become superfluid. When the two populations of 
fermions are unequal, not every particle can find a partner. Will the system neverthe- 
less stay superfluid? Here we study this intriguing question in an unequal mixture of 
strongly interacting ultracold fermionic atoms. The superfluid region vs population 
imbalance is mapped out by employing two complementary indicators: The presence 
or absence of vortices in a rotating mixture, as well as the fraction of condensed 
fermion pairs in the gas. Due to the strong interactions near a Feshbach resonance, 
the superfluid state is remarkably stable in response to population imbalance. The 
final breakdown of superfluidity marks a new quantum phase transition, the Pauli 
limit of superfluidity. 
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The study of superfluidity in an unequal mixture of two Fermi gases, or more generally 
unequal Fermi surfaces, is highly relevant to a wide area of physics. For example, this 
situation should occur in strongly degenerate quark matter in the core of a neutron star, 
where quarks of differing mass will belong to Fermi spheres of differingsize. The ground 
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and exper- 



state of such a system has been the subject of debate for decades 
imental studies are highly desirable. However, the experimental realization of unbalanced 
Fermi seas in superconductors, charged fermionic superfluids, poses extreme difficulties: To 
reach imbalanced electron densities of spin up vs spin down electrons, a natural idea would 
be to apply a magnetic field. However, magnetic fields are either fully shielded from the 
superconductor via the Meissner effect, or they enter only in the form of quantized flux 
lines or vortices. Experiments trying to study mismatched Fermi surfaces therefore have 



or 



to suppress these effects, as in experiments in heavy fermion superconductors 
quasi-two-dimensional organic superconductors In the neutral superfluid helium-3, one 
can mismatch the Fermi surfaces by a magnetic field and thus destroy inter^spin pairing. 
However, superfluidity persists due to (p-wave) pairing betwee n equal spins 
The recently discovered atomic fermionic superfluids 
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Detween equal spins |9j. 
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22j provide an exciting new possibility to explore unequal mixtures of fermions. Here, 



populations in two hyperfine states of the fermionic atom can be freely chosen. In addition, 
the (s-wave) interactions between two atoms in different states and hence the binding energy 
of atom pairs can be tuned. In equal mixtures of fermions, this tunability is being exploited 
to study the crossover from a Bose-Einstein Condensate (BEC) of molecules to a Bardeen- 
Cooper-Schrieffer (BCS) superfluid of loosely bound pairs 



zero temperature, this crossover is smooth 
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22j. At 



|25I |. the system stays superfluid even 



for arbitrarily weak interaction and no phase transition occurs. In the case of unequal 
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In the 



mixtures, the phase diagram is predicted to be much richer 
molecular limit of tight binding, all fermions in the less populated spin state will pair up 
with atoms in the other state. The resulting molecular condensate will spatially coexist 
with the remaining Fermi sea of unpaired atoms. As the repulsive interaction between 
atoms and molecules is increased, the condensate will start to expel unpaired atoms, leading 



to a phase separation of the su perfl uid from the normal phase 



picture is expected 
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to extend into the BCS-limit of weakly bound pairs. 



Here, the pairing gap A prevents unpaired atoms from entering the BCS superfluid 
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As the binding and hence the pairing gap is further reduced, A will eventually become 
small compared to the chemical potential difference 5fi = fi2 — A*i between the two spin 
states, allowing unpaired excess atoms to enter the superfluid region. Close to this point, 
superfluidity will cease to exist [3^]. In the weakly interacting BCS-limit the pairing gap is 
exponentially small compared to the Fermi energy, hence an exponentially small population 
imbalance can destroy superfluidity. 

This superfluid to normal transition is an example of a quantum phase transition, which 
occurs even at zero temperature, when all thermal fluctuations are frozen out and only 
quantum fluctuations prevail. It can also be driven by increasing the mismatch in chemical 
potentials between the two spin states to the critical value of 8/jL « A, inducing collapse into 
the normal state. In this context the phase transition is known as the Pauli or Clogston 

n 

limit of superfluidity Its exact nature, whether there is one or several first-order and/or 
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second-order transitions, however, is still the subject of current debate |6, 

In this paper we map out the superfluid region as a function of population imbalance, 
interaction strength and temperature in an ultracold fermionic gas of 6 Li atoms. The only 
direct and unambiguous signature of superfluid flow in Fermi gases so far is the presence of 
vortices j^ . By studying unequal Fermi mixtures under rotation we establish superfluidity 
for a broad range in the population mismatch. Close to the breakdown of superfluidity, 
vortices are strongly damped and difficult to observe. Therefore, we map out the full regime 
of superfluidity by determining the fraction of condensed fermion pairs in a non-rotating 



cloud 



14. 15]. 



The two experimental methods require slightly different procedures for imaging the pair 
condensate wavefunction after release from the trap. To extract the fraction of condensed 
vs uncondensed pairs, the condensate must separate well from the thermal cloud and should 
therefore remain small. For the detection of rotating clouds, the condensate should expand 
to a large size in order to magnify the vortices. In the following, we give the parameters used 
to determine the condensate fraction in parentheses after those used for vortex detection. 

In the experiment, fermionic 6 Li atoms were sympathetically cooled to degeneracy by 23 Na 
atoms in a magnetic trap 3]. The ultracold cloud was subsequently loaded into an optical 
dipole trap (waist w ~ 120 /im) at a maximum trap depth of about 8 fiK. At a magnetic bias 
field of 875 G, a variable spin-mixture of the two lowest hyperfine states (labelled |1) and 
1 2)) was created via a Landau-Zener sweep with variable sweep rate. Interactions between 
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these two states are strongly enhanced around a 300 G wide Feshbach resonance located 
at 834 G j^). At lower values of the magnetic bias field, two isolated fermions can bind 
into a stable molecule, while at higher values fermion pairs can only exist in the stabilizing 
presence of the surrounding gas. This tunability of the binding energy provides access to the 
BEC-BCS crossover physics. The spin mixture was evaporatively cooled further by lowering 
the trap depth to 1.6 fiK resulting in radial and axial trap frequencies of v r = 110 Hz and 
v a = 23 Hz, respectively. At the same time, the magnetic field was ramped to 812 G (818 
G), which is on the BEC-side of the resonance, but still in the regime of strong interaction. 
Here, l/fc^a = 0.19 (0.11), where a is the scattering length and kp is defined as the Fermi 
momentum of a non-interacting, equal spin mixture. The rather moderate evaporation still 
leaves room for thermal molecules in an equal mixture, but was chosen to ensure efficient 
cooling of highly asymmetric mixtures, avoiding spilling of large Fermi clouds. It ensured 
that the total number of atoms iV = 7 • 10 6 (N = 2.3 10 7 ) was approximately constant and 
independent of the asymmetry between the two spin states. 

For the vortex experiment, we set the spin mixture in rotation using two blue-detuned 
laser beams (wavelength 532 nm) rotated symmetrically around the cloud at angular fre- 



22| . After 800 ms of stirring, the rotating cloud was left to equilibrate 



quency f2 = 2n 70 Hz) 
for several hundred ms. 

Starting with either the rotating or the non-rotating cloud, we then varied the interaction 
strength between the two spin states in the gas by ramping the magnetic field in 100 ms (500 
ms) to several values around the Feshbach resonance (for the condensate fraction experiment, 
the trap depth was simultaneously increased to 4 /iK [v T = 192 Hz)). After 50 ms (100 ms) 
of hold time, an image of the cloud was taken following the procedure outlined in 22j ]. In 
short, after releasing the cloud from the optical trap the binding energy of fermion pairs 
was rapidly increased by ramping the magnetic field within 2 ms (200 /xs) to 690 G, in the 
far wings of the resonance on the BEC side. Here, fermion pairs were stable throughout 
further expansion. After a total of 11 ms (14 ms) of expansion (in the remaining magnetic 
saddle-point potential) an image of either state |1) or state |2) was taken. For the condensate 
fraction data, the magnetic field was suddenly switched to 800 G right before imaging. At 
this field the molecules absorb the probe light with the same strength as free atoms. The 
images revealed the center-of-mass wavefunction of the pairs and, for rotating clouds, the 
eventual presence of vortices. For the condensate fraction experiment, the 200 /xs fast ramp 
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to the BEC-side immediately after release from the trap ensured that even large condensates 
separated well from the normal, uncondensed component. Since the ramp was fast compared 
to the radial trapping period, the size of the expanded condensate was mostly governed by 
the residual mean-field interaction at 690 G, where a = 1400 ao- 

Fig. n shows profiles of the two spin states for various spin-mixtures, on the BEC- and 
on the BCS-side of the resonance. Starting with a pure Fermi sea in state |1) we see how 
gradually, for increasing numbers in the second spin state |2) first a normal (uncondensed) 
cloud of fermion pairs emerges, then a condensate peak appears within the normal cloud 
(see also Fig. EK, b). As the condensate size increases, vortices appear in the rotating cloud. 
Naturally, the largest condensates with the largest number of vortices are obtained for an 
equal mixture. Clearly, superfluidity in the strongly interacting Fermi gas is not constrained 
to a narrow region around the perfectly balanced spin-mixture, but instead superfluid flow 
is observed for large asymmetries in the populations. 

Fig. 121 summarizes our findings for rotating spin-mixtures. It displays the number of 
detected vortices vs the population imbalance between the two spin states. The vortex 
number measures qualitatively how deep the system is in the superfluid phase: The higher 
the non- superfluid fraction, the faster the condensate's rotation will damp given the non- 
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38, B9I. We therefore observe how 



vanishing anisotropy {^f-^f- ~ 1-5%) of our trap 
gradually, for decreasing interaction strength on the BCS-side, the superfluid window shrinks 
in size around the optimal situation of equal populations. 

A more detailed map of the superfluid regime as a function of interaction strength and 
also temperature (see Fig. El) was obtained from a study of condensate fractions, which were 
determined from cloud profiles such as in Fig. El Throughout the whole crossover region, 
pair condensation occurred for a broad range of population imbalances, demonstrating again 
the stability of the superfluid around the resonance. 

An intriguing property of the superfluid state with imbalanced populations is the clear 
depletion in the excess fermions of the majority component, see Fig. Et- The profiles in 
Fig. El present the axially integrated density, hence the true depletion in the 3D density 
is even stronger. The condensate seems to repel the excess fermions. This feature was 
observed after expansion at 690 G, where interactions are still strong (initially 1/kpa = 
2.0). The expansion, at least in the region around the condensate, is hydrodynamic and 



should proceed as a scaling transformation 
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421 ] . Therefore, the depletion observed 
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in ballistic expansion hints at spatial phase separation of the superfluid from the normal 
state. This effect was observed throughout the resonance region, and on resonance even 
when no magnetic field ramps were performed during expansion. However, to distinguish a 
phase separated state with equal densities in the superfluid region from more exotic states 
allowing unequal densities, a careful analysis of the three-dimensional density, reconstructed 
from the integrated optical densities, is necessary and will be the subject of a future study. 

We did not observe (by simultaneously imaging along the long and short axis) a mod- 
ulation in the condensate density as would be predicted for the FFLO state 2?| Q, Q|- 
However, this state is predicted to be favored only in a narrow region of parameter space 
and might have escaped our attention. 

The condensate fraction was determined from the minority component, which in all cases 
is very well fit by a gaussian for thermal molecules and unpaired atoms, plus a parabolic 
Thomas- Fermi-profile for the condensate. Fig. 0] shows the condensate fraction obtained 
for varying population difference and temperature, and for magnetic fields or interaction 
strengths around resonance. The data for 754 G, on the BEC-side of the resonance, 
shows condensation over almost the entire range of population imbalance. As the inter- 
action strength is increased towards resonance, the condensate fraction for equal mixtures 
grows However, for large population asymmetries it disappears. The condensation 

window shrinks further as we cross the resonance and move to the BCS-side (Fig. HJl-f). 

The temperature varied with number imbalance as indicated in the insets of Fig. 0] The 
maximum for equal mixtures at 754 G is likely due to the energy release when more molecules 
were formed. The observed critical population imbalance was only weakly dependent on 
temperature. This may reflect that the pairing gap is only a weak function of temperature, 
for temperatures well below the critical temperature for superfluidity 4^|- The critical 
imbalance at our coldest temperatures will thus essentially coincide with its value at zero 
temperature. 

On resonance, where the scattering length a diverges, the system is in the unitary 
regime 

QBE, 42, 3, 

where the only remaining energy scales of the system are the Fermi 
energies Ep i and Ep 2 of the two spin components. The breakdown of superfluidity will occur 
for a certain universal ratio of these two or equivalently, in a harmonic trap, for a certain 
critical population imbalance. We determine this universal number to be 5 C ~ ±70(3)% 
for our approximately harmonic trapping potential. This corresponds to a Fermi energy 
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difference 5E F = E F , 2 - E F)1 = ((1 + 4) 1 / 3 - (1 - SJ 1 / 3 ) e f = 0.53(3)£ F , with E F the 
Fermi energy of an equal mixture of non-interacting fermions. The standard BCS state is 
predicted 1| to break down for a critical chemical potential difference 5fi = \^2A. On res- 



onance, however, Monte-Carlo studies predict |26[ the superfluid breakdown to occur when 
5fi = 2.0(1)A = 1.0(1)E F . Only in the weakly interacting regime do the chemical potentials 
equal the Fermi energies. Quantitative agreement with the Monte-Carlo study would require 
that 8/j, « 2 5Ep. This is not unreasonable given the fact that interactions will reduce the 
chemical potential of the minority component. In a preliminary analysis, we indeed find 
close agreement with theory. 

Fig. |3] summarizes our findings. It shows the critical mismatch in Fermi energies for 
which we observed the breakdown of superfluidity as well as the pairing gap A versus 
the interaction parameter l/k F a. Far on the BEC-side of the resonance the superfluid is 
very robust with respect to population imbalance. Here, pairing is dominantly a two-body 
process: The smallest cloud of atoms in state |1) will fully pair with majority atoms in 
state |2) and condense at sufficiently low temperatures. On the BCS-side of the resonance, 
however, pairing is purely a many-body effect and depends on the density of the two Fermi 
clouds. As the density of the minority component becomes smaller, the net energy gain 
from forming a pair condensate will decrease. Even at zero temperature, this eventually 
leads to the breakdown of superfluidity and the quantum phase transition to the normal 
state. We have experimentally confirmed the qualitative picture that fermionic superfluidity 
breaks down when the difference in chemical potentials between the two species becomes 
larger than the pairing gap. 

In conclusion, we have observed superfluidity with imbalanced spin populations. 
Contrary to expectations for the weakly interacting case, superfluidity in the resonant 
region is extremely stable versus population imbalance. As the asymmetry is increased, 
we observe the quantum phase transition to the normal state, known as the Pauli limit 
of superfluidity. Our observation opens up intriguing possibilities for further studies on 
mismatched Fermi surfaces. One important aspect concerns the density distribution in 
the superfluid regime. Standard BCS theory allows only equal spin densities, which would 
entail complete phase separation of the superfluid from the normal density. More exotic 
solutions p allow superfluidity also with imbalanced densities. Equally fascinating is the 
nature of the strongly correlated normal state slightly below resonance. For sufficient 
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population imbalance we have the remarkable situation that bosonic molecules, stable even 
in isolation, do not condense at zero temperature, due to the presence of the Fermi sea. 
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FIG. 1: Superfluidity in a strongly interacting Fermi gas with imbalanced populations. The upper 
(lower) pair of rows shows clouds prepared at 812 G (853 G), where \jkpa = 0.2 (l/fepa = —0.15). 
In each pair of rows, the upper picture shows state |1), the lower one state |2). For the 812 G 
data, the population imbalance 5 = between Ni atoms in state |1) and N2 in state |2) 

was (from left to right) 100%, 90%, 80%, 62%, 28%, 18%, 10% and 0%. For the 853 G data, the 
mismatch was 100%, 74%, 58%, 48%, 32%, 16%, 7% and 0%. For different 6, the total number of 
atoms varied only within 20% around N = 7 • 10 6 , with the exception of the endpoints 5 = 100% 
(N = 1 • 10 ) and 5 = 0% (N = 1.2 • 10 7 ). The field of view of each image was 1.4 mm x 1.4 mm. 
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Population Imbalance 

FIG. 2: Vortex number vs population imbalance for different interaction strengths. Results are 
shown for 812 G or l/k F a = 0.2 (filled circles), 853 G (l/k F a = -0.15, empty circles), 874 G 
(l/fei?a = —0.3, filled triangles), 896 G (l/fc^a = —0.4, empty squares) and 917 G (1/kpa = —0.5, 
diamonds). 
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Distance from Center r [mm] Distance from Center [mm] 

FIG. 3: Radial density profiles of the two components of a strongly interacting Fermi gas mixture 
with unequal populations, a) and b) Profiles of the component in state |1) and |2) resp., originating 
from 883 G (1/fcpa = —0.27). The imaging procedure is detailed in the text. The population 
imbalance was 5 = 0% (red), 5 = 46% (blue) and 5 = 86% (green), c) Difference between the 
distributions in state |1) and |2). The total number of atoms was iV = 2.3- 10 . The clear dip in the 
blue curve caused by the pair condensate gives an indication for phase separation of the superfluid 
from the normal gas. d) Color-coded profiles of clouds prepared at three different interaction 
strengths. The profiles are azimuthal averages of the axially integrated density. The condensate is 
clearly visible as the dense central part surrounded by unpaired fermions or uncondensed molecules. 
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Population Imbalance 

FIG. 4: Condensate fraction vs population imbalance for several temperatures and interaction 
strengths. The total number of atoms N = 2.3 • 10 7 is constant to within 20% for all data-points 
(Tp = 1.9 fj,K for an equal mixture). For a given population imbalance, the uppermost curves 
for different magnetic fields are approximately isentropically connected. The different symbols 
correspond to different evaporation ramps. The average radial kinetic energy per molecule of 
thermal clouds in the minority component serves as an indicator for temperature and is shown in 
the insets for 754 G and 833 G for the coldest data. On resonance, for a population asymmetry of 
50%, we measure an energy oiks- 300 nK (circles), 345 nK (inverted triangles), 390 nK (squares), 
420 nK (triangles) and 505 nK (diamonds). The critical population imbalance 5 C for the breakdown 
of condensation at 754 G is about 5 7 c 54 ^ 96% and at 786 G it is <5j 86 as 95%. For the data at higher 
magnetic fields we determine 6 C through a threshold fit to the first three data points with non-zero 
condensate fraction for each sign of asymmetry. Although we could have used any reasonable 
threshold function, empirically, it was found that the function n c (l — |<5/<5 C | 3 ' 3 ) (n c - maximum 
condensate fraction) provided a good fit to all data points. Therefore it was used for the threshold 
fits and is shown as a guide to the eye. 
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FIG. 5: Critical difference in Fermi energies SEp between the two spin states for which the 
superfluid to normal transition is observed. 5Ep for each interaction strength and tempera- 
ture is obtained from the critical population imbalance determined in Fig. ^ using SEf/Ep — 
(1 + 5c) 1 / 3 - (1 - 5c) 1 / 3 . The symbols are defined in Fig. The line shows the expected variation 



of the pairing gap A, where the value on resonance has been taken from [26( and the exponential 
behavior in the BCS-regime, A ~ e -T/2fc F |a| wag 

assumed. While the trend of 5Ep is expected 
to follow that of A, the close agreement is coincidental. Representative density profiles illustrate 
the quantum phase transition for fixed interaction and for fixed population imbalance along the 
dashed lines. 
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